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Abstract. We study the 3D Edwards- Anderson spin glasses, by analyzing spin-spin 
correlation functions in thermalized spin configurations at low T on lattices of sizes up 
to 32 3 . We consider individual disorder samples and analyze connected clusters of very 
correlated sites: we analyze how the volume and the surface of these clusters increases 
with the lattice size. We qualify the important excitations of the system by checking 
how large they are, and we define a correlation length by measuring their gyration 
radius. We find that the clusters have a very dense interface, compatible with being 
space filling. 
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1. Introduction 

The understanding of the low temperature behavior of spin glass systems [1] is both a 
very interesting and challenging problem and one of the most relevant open questions in 
classical statistical mechanics. The mean field theory of the model has been solved [2], 
and the rigorous mathematical proof that this solution is correct [3] made the situation 
very clear. In finite dimensional spin glass however, the situation is completely open, and 
there is no general agreement on what happens in the infinite volume limit (V —> oo). 
The two main possibilities are on one side that in a finite dimension the model behaves 
as in mean field j2] (this is the replica symmetry breaking, RSB, picture), by preserving 
all or most of the very peculiar features of the mean field theory, and on the other 
side that the so called droplet picture [I] (where only two stable states are relevant for 
the critical behavior of the system) is realized. For some detailed analysis of the two 
points of view see for example [3 [6] . The TNT scenario [6] suggests an intermediate 
situation (the trivial-non-trivial picture) where the link overlap does not have a complex 
behavior; here we only estimate the usual overlap. Even if it is very plausible that the 
final answer to this crucial question will not come from numerical experiments, it is clear 
that today large scale, large volume, very accurate numerical simulations can allow to at 
least approach the relevant low T, large V, regime, and give very relevant hints about 
the physical behavior of such systems in the thermodynamic limit. The numerical 
simulations of the Janus supercomputer [HI |9j [TOj [11] (a special purpose computer [12] 
specially built to simulate very effectively spin systems with discrete variables and 
couplings) have produced a large set of well thermalized spin glass configurations on 
large (since we are considering spin glasses: the same lattice sizes would not be so large 
for, say, an Ising model with random field) systems at low T: our analysis will be based 
on these spin configurations. 

We take the approach, that has been discussed in [13], of considering local 
magnetizations and correlation functions for a given disorder sample (without averaging 
over the disorder). This approach can give a very complete picture, and allow detecting 
important features of the system. We will mainly use spin-spin correlation functions, 
k\ J j = (o-iO~j) (where the brackets indicate the thermal average, % and j are two lattice 
sites, and (J) denotes a given realization of the couplings) and build upon how much 
two sites are correlated in a given sample. 

Using the values of the kij we will build connected clusters of correlated spins. 
These clusters are related to excitations that during the dynamics tend to move in a 
coherent way: their spins are very correlated. In this sense these connected clusters 
are good proxies of droplets and studying their properties allows us to get information 
about the droplets. In ref. [H] a similar approach was taken to study excitations at 
T = 0. It is thanks to the availability of well-thermalized configurations from Janus for 
many disorder realizations that we are able to analyze the finite small T region. 
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2. Model and data analysis 

We study here the 3D Edwards Anderson spin glass with binary couplings. We analyze 
well thermalized, equilibrium spin configurations obtained with large scale numerical 
simulations by the Janus collaboration jHJ |9] based on parallel tempering [15] (see 
ref. [HI |9] for all the parameters of the numerical simulation). We deal with 0(1000) 
samples of quenched disorder and many independent spin configurations per disorder 
sample (for each lattice size), on simple cubic lattices of sizes going up to V = L 3 = 32 3 , 
for temperatures down to 0.64 T c , with periodic boundary conditions. 

For a given disorder sample and lattice site i, we form the (connected) cluster d by 
adding recursively to Ci all sites j whose correlation function ((Tj&i) is larger in absolute 
value than a given threshold T m in- When this is done, we only keep the connected 
component containing the site %. The value of the threshold Tmin is an important 
parameter of our construction: if it is very small the cluster will fill all the available 
space, and it will have no boundary, while if it is very large C\ will be small (but its 
complement, formed from all lattice sites that do not belong to Ci, will be large). We 
would expect that in the physically relevant range of values of Tmin the universal results 
of our computation (like, for example, the critical exponents) should not depend on 
Tmin-, and we will verify that this is indeed the case. 

For a given sample we will consider both the set of the V (connected) clusters Cj, 
analyzing their volume Ni and their surface Si (defined as the size of the set of sites 
which belong to the cluster and have at least one first neighboring site that does not 
belong to it), and properties obtained by averaging over the different clusters Ci (e.g. 
the cluster average volume N and average surface S). All these quantities should be 
labeled also with an index (J), denoting the disorder sample we are considering, but 
we will omit this label for the sake of a clearer notation. We will try to determine how 
quantities like TV and S scale with the lattice size. 

3. Size and surface of connected clusters 

We start by showing what happens in two representative disorder samples in figure [Ua 
and[TJb. Both figures are scatter plot that shows the cluster surface density Sj = Si/V 
versus the cluster volume density Vi = Ni/V for cluster Ci. Each point is for one of the V 
clusters one can construct when considering the disorder realization J. We show data for 
our largest lattice, with L = 32 (also, here and in the rest of this note, T = 0.7 ~ 0.64T C ). 
Different colors are for clusters reconstructed by using different values of Tmin- Scaling 
for the different values of T m in is good: data points with different colors fall on the 
same curve. The two samples we show behave quite differently: in the one on the right 
(where larger clusters are present: the maximum rij value is larger than 0.8) Sj versus 
Hi saturates and starts to decrease for rij ~ 0.7. This decrease has a clear reason, due 
to finite size effects: when clusters reach a size comparable to the lattice size (that has 
periodic boundary conditions) it fills up all the available space, and it has no space 



Correlated domains in spin glasses 4 




0.1 0.2 0.3 0.4 0.5 0.1 0.2 0.3 0.4 0.5 

N/V N/V 



Figure 1. In figurcs[I](a) and[T](b) wc plot the quantities Si/V versus Ni/V for all 
clusters Ci (i = 1, V) for two representative disorder samples. Symbols of different 
colors and shapes are used (like in figures [TJ (c) and[TJ(d)), for different values of the 
threshold T m i n : black squares are for a threshold of 0.5, cyan triangles pointing down 
0.6, blue triangles pointing up 0.7, green circles 0.8 and red squares 0.9. In figure 
[U (c) we show averages over all clusters d for a given sample (each point is for a given 
disorder realization, all disorder realizations are plotted), under the constraint A = 0.5 
(see text): in figure [T] (d) we restrict the average to percolating clusters (A = 0.5). 



to have a boundary. As — > 1 the maximum allowed surface density goes to zero. 
Because of this, we have introduced a cutoff A, and we discard from our analysis the 
clusters such that n, > A. We have also analyzed the data for different values of A 
in the range (0.25, 1.0), and our claims will turn out not to depend sensitively on the 
choice of A. In the two lower figures [Uc andHJd we show averages of the cluster volume 
and surface density (with cutoff A = 0.5) over all sites in individual realizations of the 
disorder. Each point represents one disorder sample. In the plot on the left, an average 
is done over all clusters (with A = 0.5) while on the right the average is done over 
percolating clusters only (defined as connected clusters that span the whole lattice in 
at least one direction). When one considers all clusters, the small clusters, containing 
a number of spin of order one, play an important role. Selecting and analyzing instead 
only percolating clusters one is keeping only the largest clusters (again with A = 0.5): 
in a non-critical situations, with a finite correlation length and correlations that fade 
away exponentially no such clusters would exists in the infinite volume limit. The fact 
that we do find that for increasing V a finite fraction of clusters is indeed percolating is 
a first signature of the fact that the basic, low energies domains are highly non trivial 
sets of sites. When we look at s = S/V versus n = N/V averaged over all clusters, 
we see that there is a clear curvature: this is reasonable since for small clusters S is 
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necessarily very similar to V. In the case of percolating clusters one sees far less of a 
curvature: clearly here small clusters are absent, since there is a minimal value of N 
needed for percolating. 

In order to proceed to a more quantitative assessment of the situation we need to 
go a step further. Asymptotically one expects that for large volumes N oc L dN , where 
the over-line is for the average over the quenched disorder and is the exponent that 
characterizes the asymptotic growth. We are ignoring here sub-leading corrections since 
our numerical data would not allow to fit them (we work with L = 16, 24 and 32, i.e. we 
use three data points in each fit): also it turns out that the best fits to the simple form 
with the leading scaling behavior are very good, with low values of the \ 2 - Analogously 
we define the exponent d$ by the rate of increase of the cluster surfaces with the volume, 
i.e. S<xL ds . 

By fitting our numerical data we obtain the values of the exponents that we show 
in figure [2j In the case of the ferromagnetic Ising model our procedure, as applied to the 
magnetization-magnetization connected correlation functions, would find, for T ^ T c , 
finite clusters of typical size £. We plot the values of the exponents in the two cases 
where we either consider all clusters, or only percolating clusters. A few comments are 
in order. The dependence of the exponents on the value of T m %n is very weak: they 
both stay in the range (2.8, 3.3) for values of T m in ranging from 0.5 to 0.9. In all these 
cases d^ and d$ are very similar: their difference never exceed 0.1. So, to summarize, 
both exponents turn out to be stable as a function of T m %n and A, close to the values 
3 and very close together. The low energy droplets in our system, at a value of the 
temperature deep in the broken phase, appear to have, as far as we can observe on the 
finite volume we can thermalize, a very diffused interface, compatible with being space 
filling (namely d$ — 3). 

4. Sponges 

In a situation where the system is non critical we expect connected clusters to have a 
(finite) size of the order of the correlation length: survival, in the infinite volume limit, 
of extended clusters is a signal of criticality. Such extended clusters can be defined, on 
a finite lattice, in many different ways: we will use here a very restrictive definition that 
we will denote by the word "sponge" [H] . We declare a connected cluster to be sponge- 
like if both it and at least one connected component of its complement span the whole 
lattice in the x, y and z directions. Only percolating clusters contribute to sponges. 
We use the same cutoff A = 0.5, excluding clusters that take a too large fraction of the 
lattice. This applies to C% and to clusters in the complement of Cj. 

We show in figure |3] the fraction of sponges as a function of the (linear) lattice size. 
We use different values of T m in- We compute the fraction of sponges for a given disorder 
sample and we average this fraction over the random quenched disorder. The fraction of 
sponges is sizable (of order 1/3 with the very low cutoff A = 0.5) and is never decreasing 
for increasing lattice size: the clear signal one gets from these numerical data is that 
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Figure 2. The exponents ds and d]y as a function of the correlation threshold %nn 
considering either all connected clusters (bottom) or percolating clusters only (top). 
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Figure 3. 

Fraction of sponges as a function of the (linear) lattice size. 



sponges will survive in the infinite volume limit, and the typical configurations of the 
system will be critical. It is clear that a numerical study of this kind cannot, for many 
reasons, give a clear cut answer, but in the very large volume range where the Janus 
numerical simulations allow us to study equilibrium at low T we are getting a clear hint 
towards survival of extended excitations (using different values of A would not change 
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Figure 4. £g/L versus L. The correlation length £g is denned from the non 
percolating cluster with the largest radius. 



the picture): this is what one would expect in a situation where RSB holds. The same 
situation was found by analyzing ground states of the system in ref . [H] . Here we are 
working at finite T < T c , with thermalized spin configuration. The two results are in 
some sense complementary: we are interested in a T = fixed point, so we would like 
to be at low T < T c , as close to T = as possible (in principle at T infinitesimally close 
to zero). The analysis at T = gives very interesting information, but cannot exclude 
the possibility that T = is a point with a special behavior. 

5. Typical lengths of non-extended clusters 

Last we investigate the geometrical structure of the connected clusters, trying to 
determine a typical correlation length: we summarize our results in figure HJ For 
each given disorder samples we only consider the clusters that do not percolate in any 
direction, and select the one with the largest gyration radius: 



where r k is the square distance of site k from the center of mass of the cluster, and the 
sum runs over all sites of the cluster. For fixed T and V we average £q over all disorder 
samples to obtain £g, that we plot, divided by L in figure HI We plot the results for 
different values of T m in with different symbols. 

Selecting the largest, non percolating cluster, allows us to pick up the largest 
physical length scale on a given finite lattice: the procedure is analogous to the one 
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used in a normal ferromagnet to measure a physical, finite, correlation length below the 
critical point. 

In a spin glass we have a critical temperature and all temperature values that are 
below it are possibly critical (i.e. they can have an infinite correlations length): because 
of this the analysis is more delicate than in a ferromagnet, We note first that the result 
does not depend much on T m i n and A: the emerging picture is really very stable. A 
constant value of ^q/L would imply a divergence of the correlation length typical of 
a critical point: we observe indeed a slight increase of this ratio, that we attribute to 
finite size effects. Again, all the signatures we are detecting call for a picture where a 
diverging length is determining the behavior of the system. Using the non percolating 
cluster with the largest number of sites, instead of the one with the largest radius, does 
give the same kind of behavior. 

6. Conclusions 

We have analyzed large-scale excitations of the system, estimated critical exponents, 
the probability of finding large excitations, and the gyration radius as an estimator of 
a correlation length. We have been able to work on thermalized spin configurations at 
low T on large lattices (experiments on spin glasses [16] show that we are working on 
length scales not far from the physically relevant ones). 

The emerging picture is strongly suggestive of the presence of RSB (and compatible 
with a TNT scenario): typical excitations are large and their interface is very dense, 
compatible with being, in the infinite volume limit, space filling. In order to have 
exponents compatible with the exponents we have measured, droplets should be very 
different from the excitations originally proposed in the "droplet picture" 
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